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We introduce discrete multivortex solitons in a ring of nonlinear oscillators coupled to a central site. 
Regular clusters of discrete vortices appear as a result of mode collisions and we show that their stability is 
determined by global symmetries rather than the stability of constituent vortices. Stable multivortex solitons 
support complex vortex dynamics including charge flipping and spiraling. © 2011 Optical Society of America 
OCIS codes: 190.5940, 190.6135. 



Discrete vortex solitons have received considerable at- 
tention over the past decade [1^3], realizing persistent 
currents between coupled nonlinear oscillators, such as 
optical and matter-wave waveguides. A stable discrete 
vortex requires balanced currents between each pair of 
sites in the loop. For general configurations with multiple 
vortices the balance condition should include each pair of 
neighboring vortices and thus it extends over the whole 
field, making multiple vortices globally linked [4] . Multi- 
vortex fields appear naturally in periodic media [5,6] but 
so far the only stable multivortex solitons predicted [7] 
and observed [8] in hexagonal photonic lattices consist 
of simple alternating vortex-antivortex pairs. 

In periodic media the vorticity is restricted by discrete 
symmetries [9, 10] and it is a nontrivial question if sta- 
ble discrete solitons can be constructed with asymmetric 
multiple vortices and globally linked currents. The inter- 
play between symmetry and nonlinearity is captured by 
purely discrete systems [1,2,11], such as a ring of N cou- 
pled oscillators [12], governed by the discrete nonlinear 
Schrodinger equation (DNLS). The multivortex solitons 
can be studied in this simple system if we add a central 
site [13], similar to hexagonal lattices when = 6 [7,8]. 

In this Letter, we investigate discrete vortices in a 
ring of nonlinear oscillators coupled to a central site. We 
demonstrate that multivortex linear modes and solitons 
bifurcate at the crossings of modes with different symme- 
tries, controlled by the coupling strength with the cen- 
tral site. We study in detail the heptamer with iV = 6 
and find analytically and numerically a variety of sta- 
ble multivortex solitons. By tracing the trajectories of 
individual vortices in the stable perturbed solitons we 
show the appearance of breather-like states with peri- 
odic charge flipping and spiraling of individual vortices, 
i.e. formation of helical vortex structures. 

We assume that the sites on the ring with complex 
amplitudes En are identical while the central site Eq 
can be different, characterized by the coupling strength 
C. The DNLS takes the form 

id.Eo + CY,En+ So\Eo?Eo = 0, 
id.En + En^i + En+i + CEo + Si\E„\^Er, = 0, (1) 



here and below the summation is over n = 1, N sites 
on the ring, En+N = En, and Sq, 5i = ^1 stand for 
de/focusing nonlinearities. This system conserves the to- 
tal power P =\Eo\'^ + J2 l-E'r^|^ i.e. dP/dz = 0. The ge- 
ometry for heptamer with iV = 6 is shown in Fig. 1(a). 

We look for stationary solutions with propagation con- 
stant fc, Eq = Aexp{ikz) and E„ — B„ cxp{ikz), and 
remove the global phase by setting A> 0, 

{k-doA^)A = C^i?„, 
(fc-(5i|S„|2)B„ ^ CA + B„^i+B„+i. (2) 

In the linear limit do,i = this system is solved with the 
discrete Fourier transform (DFT), i3„ = -^s^'®'", 
here Qg = 2tts/N and Eg are the Fourier amplitudes. 
First, we find vortcx-frec (real) modes with A > and 
Es=0 for ah s^N, 



k± = l± Vl + NC^ = NCFn/A. 



(3) 



The "bright" mode with k ~ k+ has all sites oscillating in 
phase, Em > 0, while the "dark" mode with k = k^ has 
the ring sites 7r-out-of-phase with the center, E^ < 0. 
The rest of the modes are constant-amplitude discrete 
vortices of charge |to| < N/2 on the ring with zero cen- 
tral site A ~ [12], i.e. with Eg^m = and the dispersion 
relation k = km — 2cosOm. Figure 1(b) shows the de- 
pendence of mode propagation constants on the coupling 
C for a heptamer. 

At specific values of the coupling C the k- curve inter- 
sects with k,n lines in Fig. 1(b), allowing the dark mode 
to coexist with vortex modes in a stationary superposi- 
tion. The collision of two modes creates novel states with 
distinct topologies, or multivortex modes, with multiple 
vortices residing in different elementary triangular cells, 
as we show schematically in Figs. 1(q;,/3,7). While the 
discrete vortices without the central site, C = 0, obey 
the symmetry En+i — Enexp(iQrn), this symmetry is 
broken for multivortex modes. As a result, m no longer 
uniquely defines the topology of the mode. For example, 
in Fig. 1(b) the intersection of two modes fc_ and k2 sup- 
ports two different configurations, (a) with to = 2 and 
(/3) with TO = 0. 



1 




k 









































Fig. 1. (Color online) (a) Heptamer configuration of non- 
linear oscillators, (b) Wavenumbers of linear modes vs. 
central coupling C. (a, /3, 7) Multivortex modes at cross- 
ings shown by dots in (b); the arrows show non-zero 
currents and the numbers (—1,0,-1-1) indicate topologi- 
cal charges m„ in each triangular cell. The total charge 
TO = ^ m„ is TO = 2 for mode {a) and m = for modes 
(A 7)- 

We now examine the nonlinear case, starting with con- 
stant amplitude solitons, B„ = B cx'p{iip„) and B > 0. 
With this Ansatz the exact solutions can be obtained by 
applying DFT, exp(i^/>„) = exp(i6sn), thus 

Bbsik-6iB^ -2 cosSs) ^0, s^N 

BbNik- SiB'^ -2) = CA, 

BbNNC ={k- doA^)A. (4) 

As in the linear case above, there are vortex-free solitons 
with bs^N = 0, as well as the discrete vortex solitons 
with A = and fc = fc^ = 2 cos 9,„ -I- SiB^ [12]. 

The discrete multivortex solitons appear, similar to 
the linear multivortex modes, as superpositions of a real 
mode with amplitude 67V and two vortex modes. 



J4', 



(5) 



with k = kjn and complex amplitudes b±jn constrained 
by |exp(iV;„)| = 1. We find such solutions for A^/|m| = 
2,3,4 only. In particular, the N = 2\m\ solitons with 
alternating charges, as in Fig. 1(7), have been observed 
experimentally in hexagonal photonic lattices [8] . 

Solutions of particular interest in application to the 
heptamer in Fig. 1 are the ones with N ~ 3|m|, realizing 
topological configurations in Figs. l(a, (3). Periodicity of 
Eq. (5) requires V'j+a ~ ''Pj with only three indepen- 
dent phases, j — 1,2,3. Thus, instead of Eq. (5) with 
two complex parameters b±m, it is convenient to use the 
identity J2 = Nb^- We derive from Eqs. (4) 



k = SiB^ - 1 = 6aA^ 



-CA/B, 



(6) 



The real and imaginary parts of Eq. (6) form two equa- 
tions for three variables V'i,2,3- Therefore, one of three 
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Fig. 2. (Color online) Existence and stability diagrams 
for solitons with (5o,i ~ 1 and C — 0.6 (a), C = 1 (b), 
and C = 1.3 (c). Solid/dashed lines indicate linearly sta- 
ble/unstable solitons. Blue, black, and green curves cor- 
respond to the dark, vortex, and multivortex solitons. 
Columns (d,c) show trajectories of vortices with charges 
±1 (red and blue lines) for stable perturbed multivortex 
solitons (a,/?, 7) indicated by dots in (a,b). Initial per- 
turbations of solitons are Eq = 0.6^ for (a), Eq = 0.8^ 
for (/3), and \En\ = 0.8J5 for (7). 



phases is a free parameter, leaving the total power and 
the Hamiltonian invariant. The total charge and posi- 
tions of vortices within the hexagon in Figs. l(a,/3) de- 
pend on this free parameter as they trace a circle for 
CA/B < 1 with net charge ±2, while for CA/B > 1 they 
trace a hexagon with zero net charge. We demonstrate in 
Fig. 2 that such degeneracy of stationary solutions leads 
to their nontrivial dynamics with spiraling vortices. 

In Fig. 2(a-c) we present existence and stability dia- 
grams for constant-amplitude solitons in the heptamer. 
At small coupling, such as C = 0.6 in Fig. 2(a), the 
multivortex solitons (a, 7) bifurcate from vortices with 
|to| = 2,3, respectively. As the coupling is gradually in- 
creased, the collision of the linear modes marked (a, /3) in 
Fig. 1(b) occurs, and the dark, vortex, and multivortex 
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Fig. 3. (Color online) Examples of stable (solid lines) and 
unstable (dashed lines) asymmetric multivortex solitons 
for = 6, (5o.i = 1, and C = 1. Each color corresponds 
to a distinct vortex configuration shown in the insets. 

solitons bifurcate (P > 0) from this point (not shown). 
For stronger coupling, such as C = 1 in Fig. 2(b), the 
multivortex soliton (/?) with m — Q bifurcates from the 
dark soliton with With further increase to C = 1.3 
in Fig. 2(c) the bifurcation for (7) soliton changes the 
"parent" branch from k2 in (a,b) to /c_ in (c). 

We test the stability of three different solitons, marked 
by dots in Fig. 2(a,b), by numerical simulations of ini- 
tially perturbed stationary solutions. The results are pre- 
sented in Fig. 2(d,e) as the trajectories of vortices, recov- 
ered from numerical solutions to Eq. (1) as described in 
Rcf. [12]. At small propagation distances in the column 
Fig. 2(d) we observe small-scale "vibrations" of vortices, 
with the frequencies given by the eigenvalues of internal 
modes, as expected for linearly stable solitons. However, 
at large propagation scales in the column Fig. 2(e) we 
observe formation of adiabatically slow superstructures, 
in particular with vortices spiraling in (a, /3) and mov- 
ing between triangular cells, thus realizing cyclic charge- 
flipping. In contrast, the vortices in the excited multivor- 
tex soliton (7) oscillate in radial direction and flip their 
charges simultaneously, via an avoided crossing of vortex 
lines, similar to the vortex breathers in Ref. [12]. 

We found that the parent branch always loses stability 
at the multivortex bifurcation point. Thus the appear- 
ance of multiple regions of stability and instability in 
Fig. 2(a-c) suggests the existence of stable asymmetric 
multivortex solitons. Several examples of such solitons, 
preserving one axis of symmetry, are obtained numeri- 
cally from Eqs. (2) and presented in Fig. 3. The stability 
of these stationary states is determined by their global 
symmetries [14], for example the soliton with a single 
vortex within hexagon [Fig. 3, black lines] has a stability 
window, while its symmetric counterpart, with the vor- 
tex located on the central site, is always unstable [12]. 
We also find completely asymmetric stationary states, 
but in general they appear to be unstable. 

In conclusion, we have introduced multivortex solitons 
in an oligomer cluster of coupled nonlinear oscillators. 
Such linear and nonlinear states appear at the intersec- 
tion of vortex and vortex-free modes. The perturbation 



of stable multivortex solitons leads to the excitation of 
breather-like states with nontrivial vortex dynamics, in- 
cluding periodic charge flipping and spiraling. 

We note that the system considered here resembles 
the heptamer clusters of plasmonic nanoparticles [15- 
17]. Plasmonic nanoclusters exhibit sharp Fano res- 
onances [18] in their transmission spectrum, as the 
consequence of interference between dark and bright 
modes [15] introduced [16] and tuned [17] by the cou- 
pling C with the central site. Our results suggest that 
the collision of plasmonic modes can lead to the excita- 
tion of multiple persistent currents in metallic nanoclus- 
ters, with additional Fano-like resonances corresponding 
to different multi-vortex configurations. 

This work is supported by the Australian Research 
Council. 
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